Kinetics of the collisionless expansion of spherical nanoplasmas 
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y^ ' the behavior of the ion energy spectrum, which is monotonic only for high electron temperatures, 

^ . otherwise exhibiting a local peak far from the cutoff energy. 



Abstract 

The collisionless expansion of spherical plasmas composed of cold ions and hot electrons is ana- 
lyzed using a novel kinetic model, with special emphasis on the influence of the electron dynamics. 
Simple, general laws are found, relating the relevant expansion features to the initial conditions 
of the plasma, determined from a single dimensionless parameter. A transition is identified in 
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Recent experiments on the interaction of ultraintense laser pulses with atomic and molec- 
ular clusters [l|, |^ have shown the possibility of accelerating ions to energies of interest for 
many applications, such as nuclear fusion and X-ray generation. The interpretation of these 
experiments and the need of controlling critical features of the ion acceleration (in particular, 
their energy spectrum) require a theoretical insight of the laws governing the collisionless 
expansion of a spherical plasma driven by hot electrons. A detailed knowledge of the kinet- 
ics of the expansion is also necessary for particular applications, such as the biomolecular 
imaging with ultrashort X-ray pulses 3|, where sample damage before the imaging time 
must be avoided. Moreover, the strong interplay between ions and electrons when using 
high-intensity lasers can be used to control the expansion of large clusters, by tailoring the 
ion phase space jjj. At present, analytical solutions for a spherical expansion exist only 
for ideal cases, such as the Coulomb explosion (CE) p\ of a pure ion plasma, which occurs 
when all the electrons are suddenly swept away from the cluster by the laser field. In op- 
posite conditions, hydrodynamic models have been proposed p to estimate basic features 
of expansions far from the CE regime, while a kinetic solution for the adiabatic expansion 
of plasma bunches into a vacuum has been derived in the quasineutral limit 7]. However, 
such extreme situations are hardly met in experiments, in particular for large clusters, as 
the laser intensity is insufficient to drive a pure CE but high enough to produce relevant 
charge buildups within the clusters: in such cases, the expansion is described correctly only 
by kinetic models, based on the Vlasov-Poisson equations. In this Letter, the long-term 
dynamics of the expansion of spherical plasmas is analyzed by using a novel Lagrangian 
model, which allows a self-consistent, kinetic description of the radial motion of the ions 
and of the three-dimensional motion of nonrelativistic electrons, with high accuracy and low 
computational effort. The important case of initially-Maxwellian electrons is investigated in 
detail, with special emphasis on the influence of the initial conditions. Simple relationships 
are deduced for the most important physical quantities, valid for a wide range of character- 
istic parameters of the expansion. These empirical formulae are surprisingly accurate and 
suggest the possibility of a rigorous analytical derivation. Furthermore, the study of the 
asymptotic phase of the expansion shows that the behavior of the ion energy spectrum for 
realistic conditions differs greatly from that of a pure CE. Finally, the importance of these 
results to interpret experimental data is discussed. 

For the present analysis, the electrons are assumed to be heated instantaneously, as upon 



irradiation with an ultrashort laser pulse, and their initial distribution is assumed to be 
Maxwellian, with temperature Tq. The expansion process is divided in two stages: a rapid 
expansion of the electrons, which leads to an equilibrium configuration before the ions move 
appreciably, and a subsequent, slower expansion of the plasma bulk. The analysis of the 
early expansion (characterized by density oscillations that are damped on a time scale much 
faster than the ion motion) is necessary in order to evaluate the equilibrium distribution of 
the electrons, which is used as initial condition for the bulk expansion. The dynamics of 
the plasma is analyzed self-consistently by following the motion of the ions, along with the 
evolution of the electron distribution, which is described as a sequence of ergodic equilibrium 
configurations, so that the only independent variable related to the electrons is their total 
energy e = |meV^ — 6$ (the electrostatic potential $ is set to zero at infinity, so that e < 
for trapped electrons). A Lagrangian approach is used for both ions (which move radially, 
starting from the initial position rg, with zero velocity) and electrons (whose energy e evolves 
in time starting from the initial value eo). The model determines the ion trajectories ri(ro, t), 
the electron energies e(eo,t), the ion density ni(r, t), the electron density ne{r,t), and the 
potential $(r, t), according to the set of equations 
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where m^ is the ion mass, and Z the ion charge state. The expansion dynamics is determined 
once the initial ion density Ui^q and the electron energy distribution pe,o are given. In Eq. 
()ld|) . the electron density is expressed as the sum of the number of electrons pe(e)de = 
Pe,o(eo)deo having energy in [e, e + de], multiplied by the probability for an electron with 
energy e to be found at the radius r, according to the ergodic distribution. Due to the large 
mass disparity between ions and electrons, the ergodic invariant [s| for the electrons [i.e., 
the phase-space volume enclosed by the surface of equation ^meV^ — e$ = e, as defined 
in Eq. (|Te))] is conserved; this determines the electron energy distribution at any time 



during the expansion of the plasma bulk. Even though Eq. (|Tc|) is written under the 
hypothesis of no ion overtaking {dri/dro 7^ 0) p, the model can be easily generalized to 
include many-branched shock shells j^ and different ion species. As can be noticed by 
writing Eqs. (1) in dimensionless form, the dynamics of the expansion depends on the single 
parameter Tq = Z/cbTo/cce = SA^Z-Rg, being Ad the Debye length for the electrons, and 
ecE = ZeQo/Ro the maximum ion energy attainable from the CE of a uniformly-charged 
sphere of ions, with radius Rq and total charge Qq. 

The study of the bulk expansion does not require the detailed knowledge of the early 
expansion transient, but only its final equilibrium. Equations (jldjl and (fTe|l are valid only 
for sufficiently-smooth variations of $(r, t) in time, a condition which is not met in the early 
stage when the electrons are suddenly allowed to expand (as if a rigid wall, initially confining 
the hot electrons, were instantaneously brought to infinity). However, these equations can 
still be used by introducing a rigid potential barrier and gradually moving it outward, 
provided that special care is taken to avoid any energy exchange between the electrons and 
the expanding wall [this is done with a suitable adjustment of Eq. f|Te|) ]. which would lead 
to an overestimate of the electron cooling. The validity of this procedure has been tested 
by using reference results from particle-in-cell simulations [9J. In the comparisons, a cluster 
of radius Rq = 32 nm was considered, with electron density 2.3 x 10^^ cm~'^, for different 
electron temperatures spanning the range 1-10 keV. In Figure ^ the equilibrium electron 
density and the electric field are compared for Tq = 1 keV (Tq = 7.2 x 10^^) and 10 keV 
(To = 7.2 X 10^^), showing the excellent agreement between the different calculations. During 
the fixed-ion expansion, 5% of the electrons leave the ion core for Tq = 7.2 x 10^'^, while 
the percentage rises to 38% when Tq = 7.2 x 10"^. Figure |21 shows the equilibrium charge 
buildup within the ion core, AQ, and the total kinetic energy of the electrons, S (accounting 
for all confined electrons, inside and outside the ion core), as a function of Tq. Considering 
that the analytical results [li| for semi-infinite planar expansions must be recovered in the 
spherical case for Tq <^ 1 (as Ad <^ -Ro), a simple fit for AQ has been found, in the form 

Ag/go = -^2.6o(v^Tj^'), (1) 

where J^^{x) = s/ (1 + x^) . As shown in Fig. |2l the accuracy of the fit is excellent even 
when more than 90% of the electrons leave the ion core. A similar fit holds for £, as 

^/^o= 1-^3.35(1.86 To'/'), (2) 



where £q is the initial thermal energy of the electrons. In the low-temperature limit (Tq *^ 1), 
5qs. (IH) and reduce to AQ/Qq ~ a/O/cTo (as in the semi-infinite planar expansion 
lOJ l) and £/£o — 1 — 1-86 Tg , respectively. 

Once the initial equilibrium distribution of electrons is determined, thus setting the initial 
conditions of the plasma, the expansion can be analyzed by solving Eqs. (1). The early phase 
of the bulk expansion is responsible for the main part of the ion acceleration. In this stage, 
the inner ions expand much slower than the outer ions: for Tq <^ 0.1, the plasma core stays 
initially still, while a rarefaction front propagates inward until it finally encompasses the 
whole ion distribution. On the contrary, for Tq > 0.1, all the ions are promptly involved in 
the expansion. The different behaviors are depicted in Fig. El where the radial trajectories 
of isodensity points are plotted, for the same cases of Fig. ^ The case of Fig. Eb represents 
a reference situation of practical interest, in which the clusters undergo a rapid expansion 
in a hybrid regime, far from both the hydrodynamic limit and the CE regime. In such 
situation, previous models do not provide an accurate description of the expansion, because 
a self-consistent study of the electron and ion dynamics is required. As illustrated in Fig. 
m (for the scenario of Fig. Ob), the electrons rapidly cool down and the charge buildup 
within the ion front decreases, until a ballistic regime is reached for both species IHi]. In 
general, the electron dynamics strongly affects the energy spectrum of the ions, by reducing 
the repulsive electric field. Consequently, the final energies of the ions are lower with respect 
to the CE case. In fact, the asymptotic energy e^o of an ion starting at tq is given by 

eoo(ro) _ g(ro,0) ^ f^ 1 dq{n{ro,t),t) ^^ 



Ze ro Jo ri(ro,t) dt 

where g(r, t) is the net charge buildup enveloped by a sphere of radius r at time t. The 
integral term (vanishing for a CE) accounts for the energy loss due to the decrease of the 
positive charge buildup experienced by each ion along its trajectory. As expected, the shape 
of the asymptotic energy spectrum of the ions and its cutoff energy depend strongly on the 
initial conditions, as shown in Fig. El for Tq = 7.2 x 10^^, the spectrum exhibits a maximum 
at 12% of the energy cutoff, 0.27ecE- Such distribution is qualitatively and quantitatively 
different from the asymptotic spectrum of a pure CE, |(e/ecE)^''^, with cutoff energy ecE- 
As Tq increases, the distribution flattens until the maximum disappears for To = 0.5 (a 
condition close to the CE regime, with cutoff above 0.7ecE), which can hence be taken as a 
lower bound for the validity of the CE model. The cutoff ion energy emax as a function of Tq 



is shown in Fig. along with the energy corresponding to the peak in the spectrum, epeak- 
The cutoff energy admits a simple fit in the form 



^1.43 (2.28 To'/') ecE, (4) 



which, for Tq ^ 1, reduces to e^ax — 2.28 Tq ecE, while, for To < 0.5, epeak exhibits a simple 
power-law behavior, as ep^ak = 0.3TQ'^ecE- 

In recent experiments J2] on the interaction of intense lasers (10^^ — 10^^ W/cm^) with large 
hydrogen clusters (up to 2 x 10^ atoms), the measured ion spectra exhibit a local maximum, 
which has been explained, in the framework of a CE model, as an effect of the distribution 
of cluster sizes. The present analysis shows that nonmonotonic energy spectra can also 
arise from single-cluster effects, indicating that a maximum can appear even for narrow 
distributions of cluster radii (for example, with large clusters containing 10^ — 10^ atoms). 
In such cases, the spectrum features (specifically, the laws for epeak and emax presented here) 
also provide an estimate for the initial electron temperature. 

In conclusion, the dynamics of the electron-driven expansion of spherical nanoplasmas 
has been fully analyzed using a new kinetic model. Simple, accurate laws have been derived 
for the general properties of the expansion, which are valid for any value of initial electron 
temperature (as long as relativistic effects are negligible). The study also revealed peculiar 
features of the ion energy spectrum (namely, the presence of a local maximum and the 
transition from nonmonotonic to monotonic behavior) that are not taken into account by 
simplified models and that can be important for the interpretation of experiments, where 
single-cluster effects are relevant. Moreover, the strong dependence on the initial conditions 
indicates the possibility of tailoring the ion spectrum, thus improving control over the plasma 
expansion. 
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FIG. 1: Electron density (black) and electric field (gray) for (a) Tq = 7.2 x 10~^ and (b) Tq = 
7.2 X 10^^, after the early electron expansion with immobile ions. Lines refer to the present theory 
and markers to particle- in-cell simulations. 




Initial electron temperature [eQ /R j 



FIG. 2: Charge buildup AQ (gray) and total electron kinetic energy £ (black) as functions of Tq, 
after the early electron expansion with immobile ions. Circles refer to the present theory, while 
solid lines refer to the fit laws of Eqs. (^ and @. Dashed lines show the power-law behavior of 
AQ and £ in the low-temperature regime. 
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FIG. 3: Radial trajectories of the isodensity points corresponding to 98%, 75%, 50%, and 25% of 
the initial value (curves 1 through 4), for (a) Tq = 7.2 x 10^^ and (b) Tq = 7.2 x 10^^. Times are 
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FIG. 4: Evolution of the electronic charge within the ion front (gray), and of the total kinetic 
energy of the electrons (black), for Tq = 7.2 x 10~^. 




0.2 0.4 0.6 
Ion energy [ZeQ /R j 



FIG. 5: Asymptotic ion energy spectra (black) for Tq = 7.2 x 10 ^, 0.2, and 0.5. The gray curve 
refers to the pure CE case. 
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FIG. 6: Cutoff ion energy (black) and position of the maximum in the ion energy spectrum (gray) 
as a function of Tq: circles refer to the present theory, while the solid lines refer to the fit laws in 
the text. The dashed line represents the power-law behavior of emax for Tq <^ 1. 
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